BAB4R-CONF-04/039 
SLAC-PUB-10655 



o 
o 



Measurement of the Branching Fractions and CP Asymmetries of 

o 

(CP)- 



=5 
< 



; B~ — ► D? cp) K~ Decays with the BABAR Detector 



The BABAR Collaboration 

> 

00 ■ February 7, 2008 

O 
OO 
O 

"-»«. ' Abstract 
X 

We present a study of B~ — > D® Cp ^K~ decays, where D® cp ^ is reconstructed in flavor (K~ir + ), 
CP-even (K~ K + ,ir~7r + ) and CP-odd (K®tt ) eigenstates, based on a sample of about 214 million 
T(AS) — > BB decays collected with the BABAR detector at the PEP-II e + e _ storage ring. Along 
? Cabibbo-suppressed B~ — > D® Cp ^K~ decays we reconstruct also the Cabibbo-favored 
B~ — ► D®Qp^ir~ decays. We measure the double ratio of branching fractions 

B(B-^D° CP+ K-)/B(B-^D° CP+ 7T-) 



R- = 

and the CP asymmetries 

A-CP+ = 



B(B-->D°K-)/B(B-^>D°ir-) 
0.87 ± 0.14(stat) ± 0.06(syst), 



B(B-^D° cp _K-)/B(B-^D° cp _it-) 

B(B-^D°K-)/B{B-^D°ir-) 
0.80 ± 0.14(stat) ± 0.08(syst), 



B(B-^D° CP+ K-) - B(B+^D° CP+ K+) 
B(B-^D% + K-)+B{B+^Dl P+ K+) 
0.40 ± 0.15(stat) ± 0.08(syst) 



A C p- = 



B(B-^D° CP _K~) - B(B+^D° CP _K + ) 
B{B-^Dlp_K-) + B{B+^D Q CP _K+) 
0.21 ± 0.17(stat) ± 0.07(syst). 



All results are preliminary. 
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1 INTRODUCTION 



A theoretically clean measurement of the angle 7 = arg(— V u dV* b /V c dV* b ) can be obtained from the 
study of B~— >.d(*) _K"(*)~ decays by exploiting the interference between the b — > cms and 6 — > ucs 
decay amplitudes The method originally proposed by Gronau, Wyler and London is based 
on the interference between B~ — > D°K~ and B~ — > D°K~ when the D° and P>° decay to CP 
eigenstates. 

We define the ratios P and Rcp± of Cabibbo-suppressed to Cabibbo-favored branching fractions 

B(B-^ cp±) 7r-)+B(B+^ J D° Cp±) 7r+) 1 j 

with the neutral D meson reconstructed in non-CP (D°) or CP-even/odd eigenstates (D^p^_) 
channels, and the direct CP asymmetry 

= B(B-^D° CP± K~) - B{B+^D% ± K+) 
CP± B(B-^D° CP± K~)+B(B+^D° CP± K+) " U 

Neglecting the P>° - P>° mixing and the ratio 7V = A(B — » D°tt )/A(B — > P>°7r ) of the 
amplitudes of the P~ — > D°ir~ and P~ — > D°-k~ processes (\r^\ < 0.02), we can write R± = 
Rcp±/R = 1 + r 2 ± 2rcos<5cos7 and Aqp± = ±2r sin 5 sin 7/(1 + r 2 ± 2r cos <5 cos 7) . Here r = 
\A(B~ — > P>°X~)/A(P _ — > P>°i^ _ )| is the magnitude of the ratio of the amplitudes for the 
processes B~ — > D°K~ and P _ — > D°K~, expected from theory to be about 0.1 - 0.2, and 5 is 
the relative strong phase between these two amplitudes [T]. The measurement of R± and Aqp± 
allows one to constrain the three unknowns r, 5 and the CKM angle 7. In this paper we present 
the measurement of R± and Acp±- 



2 THE BABAR DETECTOR AND DATASET 

The measurements reported in this paper have been obtained from a sample of about 214 million 
T(45*) decays to BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy 
B factory. The BABAR detector is described in detail elsewhere j2]. Charged-particle tracking 
is provided by a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). For 
charged-particle identification, ionization energy loss in the DCH and SVT, and Cherenkov radia- 
tion detected in a ring-imaging device (DIRC) are used. Photons are identified by the electromag- 
netic calorimeter (EMC), which comprises 6580 thallium-doped Csl crystals. These systems are 
mounted inside a 1.5-T solenoidal superconducting magnet. The segmented flux return, including 
endcaps, is instrumented with resistive plate chambers (IFR) for muon and K® identification. We 
use the GEANT 3 j software to simulate interactions of particles traversing the detector, taking 
into account the varying accelerator and detector conditions. 



3 ANALYSIS METHOD 

We reconstruct B~ -^D°h~ decays, where the prompt track h~ is a kaon or a pion. Reference to the 
charge-conjugate state is implied here and throughout the text unless otherwise stated. Candidates 
for D° are reconstructed in the CP-even eigenstates 7r~7r + and K~K + , in the CP-odd eigenstate 
Kgir , and in the non-CP flavor eigenstate K~n + . K® candidates are selected in the 7r~7r + channel. 
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The prompt particle h~ is required to have momentum greater than 1.4 GeV/c. Particle ID 
information from the drift chamber and, when available, from the DIRC must be consistent with 
the kaon hypothesis for the K meson candidate in all D° modes and with the pion hypothesis for 
the tt^ meson candidates in the _D ^7r~7r + mode. For the prompt track to be identified as a pion or 
a kaon, we require that at least five Cherenkov photons are detected to insure a good measurement 
of the Cherenkov angle. We reject a candidate track if its Cherenkov angle is not within 3a of the 
expected value for either the kaon or pion mass hypothesis. We also reject candidate tracks that 
are identified as electrons by the DCH and the EMC or as muons by the DCH and the IFR. 

Photon candidates are clusters in the EMC that are not matched to any charged track, have a 
raw energy greater than 30 MeV and lateral shower shape consistent with the expected pattern of 
energy deposit from an electromagnetic shower. Photon pairs with invariant mass within the range 
115-150 MeV/c 2 (~3cr) and total energy greater than 200 MeV are considered ir° candidates. 
To improve the momentum resolution, the 7r° candidates are kinematically fit with their mass 
constrained to the nominal tt° mass [I]. 

Neutral kaons are reconstructed from pairs of oppositely charged tracks with the invariant mass 
within 10 MeV (~3<r) from the nominal K° mass. We also require that the ratio between the flight 
length distance in the plane transverse to the beams direction and its uncertainty is greater than 



The invariant mass of a D° candidate, m(D 1 ), must be within 3a of the D° mass. The D mass 
resolution a is about 7.5 MeV in the K~ir + , K~K + and ir~ir + modes, and about 21 MeV in the 
Kgir mode. Selected D° candidates are fitted with a constraint to the nominal D° mass. 

We reconstruct B meson candidates by combining a D° candidate with a track h~. For 
the K~7r + mode, the charge of the track h~ must match that of the kaon from the D° me- 
son decay. We select B meson candidates by using the beam-energy-substituted mass m-Es = 



J (E* 2 /2 + Pi • pb) 2 /E 2 — p 2 B and the energy difference AE = E* B — E*/2, where the subscripts i 

and B refer to the initial e + e~ system and the B candidate respectively, and the asterisk denotes 
the center-of-mass (CM) (Y(AS)) frame. The tties distributions for B~ ^D°h~ signal events are 
Gaussian distributions centered at the B mass with a resolution of 2.6 MeV/ c 2 , which does not 
depend on the decay mode or on the nature of the prompt track. In contrast, the AE distributions 
depend on the mass assigned to the prompt track. We evaluate AE with the kaon mass hypothesis 
so that the distributions are centered near zero for B~—>D°K~ events and shifted, on average, 
by approximately 50 MeV for B~^D°tt~ events. The AE resolution depends on the momentum 
resolution for the D° meson and the prompt track h~, and is typically 17 MeV for all the D° decay 
modes. All B candidates are selected with tbes within 2.5a of the mean value and with AE in the 
range -0.15 < AE < 0.18 GeV. 

To reduce background from continuum production of light quarks, we construct a Fisher dis- 
criminant based on the following quantities: (i) the scalar sum of the momenta of all charged and 
neutral particles (exluding the B decay products) flowing into nine concentric cones centered on 
the B candidate thrust axis in the CM frame; (ii) the normalized second Fox- Wolfram moment [3], 
i?2 = H2/H0, where Hi is the /-order Fox- Wolfram moment of all the charged tracks and neutral 
clusters in the event; (iii) | cos 9t\, where 0t is the angle between the thrust axes of the B candidate 
and of the remaining charged tracks and neutral clusters, evaluated in the CM frame; (iv) | cos 9g\ , 
where 6b is the polar angle of the B candidate in the CM frame; (v) | cos9h e i(D°)\, where 6h e i(D°) 
is the angle between the direction of one of the decay products of the D° and the direction of flight 
of the B, in the D° rest frame. Each cone in (i) subtends an angle of 10° in the CM and is folded 
to combine the forward and the backward intervals. A cut on the value of the Fisher discriminant 
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rejects more than 90% of the continuum background while retaining 77% of the signal in the K~ir+, 
K~ K + and K^ir modes and 65% in the it~it + channel. 

Multiple >.D /i~candidates are found in about 4% of the events for the K^tt and in less 
than 1% of the events for the other D° decays. In these events a x 2 is constructed from m(7r°) (for 
KgTT only), m(D°), and m^s and only the candidate with the smallest x 2 is retained. The total 
reconstruction efficiencies, based on simulated signal events, are about 33%(K~tt + ), 28%(K~ K + ), 



26%(7r-7r+) and l7%(K° s ir°). 

The main contributions to the BB background come from the processes B^D*h (h = tt,K), 
B-->D°p- and mis-reconstructed B~^D°h~. For D° -> K~K+, D° -> tt-tt + and D° -> K° s ir° 



be considered, since they have the same AE and mES distribution as the D K~ signal. Their 
yields are estimated from the existing measurements jlj E] and subtracted from the B~^D°K~ 
signal yields. 

For each D° decay mode an extended unbinned maximum-likelihood fit to the selected data 
events determines the signal and background yields rtj (i = 1 to M, where M is the total number 
of signal and background channels). Two kinds of signal events, B~^D°ir~ and B~^D°K~ , are 
considered, and four kinds of backgrounds: candidates selected either from continuum or from BB 
events, in which the prompt track is either a pion or a kaon. 

The input variables to the fit are AE and a particle identification probability for the prompt 
track based on the Cherenkov angle 6c, the momentum p and the polar angle 9 of the track. The 
extended likelihood function C is defined as 



where ./V is the total number of observed events. The M functions Vi(AE, 6c', at) are the probability 
density functions (PDFs) for the variables AE,9c, given the set of parameters <3j. Since these 
two quantities are sufficiently uncorrelated, their probability density functions are evaluated as a 
product V t = Vi{AE; a { ) x Vi{6 c ; a t ). 

The AE distribution for B~ -^D°K~ signal events is parametrized with a Gaussian function. 
The AE distribution for B~^D°tt~ is parametrized with the same Gaussian used for B~^>D°K~ 
with a relative shift of the mean, computed event by event as a function of the prompt track 
momentum, arising from the wrong mass assignment to the prompt track. The offset and width of 
the Gaussian are kept floating in the fit and are determined from data together with the yields. 

The AE distribution for the continuum background is parametrized with a linear function 
whose slope is determined from off-resonance data. The AE distribution for the BB background is 
empirically parametrized with the sum of a Gaussian and an exponential function when the prompt 
track is a pion, and with an exponential function when the prompt track is a kaon. The parameters 
are determined from simulated events. 

The particle identification PDF is obtained from a pure control sample of kaons and pions 
produced in the decay chain D* + — > D°ir + (D°^K~tt + ), selected using kinematical information 
only, without any inputs from the BABAR particle identification system. The parametrization of the 
particle identification PDF is performed by fitting with a Gaussian distribution the background- 
subtracted distribution of the difference between the reconstructed and expected Cherenkov angles 
of the selected kaons and pions. 



decays, the peaking backgrounds B — >K K + K , B — >K tt + tt and B 



Kq7t°K must also 
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4 PHYSICS RESULTS AND SYSTEMATIC STUDIES 



The results of the fit are summarized in Table ^ Figure ^ shows the distributions of AE for 
the K~it + , CP± and CP— modes after enhancing the B —* D°K purity by requiring that the 
prompt track be consistent with the kaon hypothesis. This requirement is about 95% efficient for 
the B~—*-D°K~ signal while retaining only 4% of the B~^-D°tt~ candidates. The projection of a 
likelihood fit, modified to take into account the tighter selection criteria, is overlaid in the figure. 

Table 1: Results of the B~^D°K- and B'^D ^ yields from the maximum-likelihood fit on 
data. 



D° mode 


N(B -» D°tt) 


N(B -> D°K) 


N{B~ -> D°K~) 


N{B+ -» 


K-7T+ 


11930 ± 120 


897 ± 34 


441 ± 24 


456 ± 25 


K-K + 


1093 ± 37 


75+ 13 
'°-12 


54^ u 


22^ 




345 ± 22 


18 ±7 


12 ±5 


7+5 
' -4 




1248 ± 40 


76±ig 


46^ u 


30+g 



The double ratios R± are computed by scaling the ratios of the numbers of B~^>D°K~ and 
B~ ^D°tt~ mesons by correction factors (ranging from 0.997 to 1.020 depending on the D° mode) 
that account for small differences in the efficiency between the B~^-D°K~ and >D°7r _ selec- 
tions, estimated with simulated signal samples. The results are listed in Tabic [2j 

The direct CP asymmetries Acp± for the B^ — > D^p.K^ decays are calculated from the 
measured yields of positive and negative charged meson decays and the results are reported in 
Table |H 

Table 2: Measured double branching fraction ratios R± and CP asymmetries Aqp± for different 
D° decay modes. The first error is statistical, the second is systematic. 



D° decay mode 


Rcp/R 


A C p 


K~K+ 


0.92 ±0.16 ±0.07 


0.43 ±0.16 ±0.09 




0.70 ± 0.29 ± 0.09 


0.27 ±0.40 ±0.09 


CP-even combined 


0.87 ±0.14 ±0.06 


0.40 ±0.15 ±0.08 




0.80 ±0.14 ±0.08 


0.21 ±0.17 ±0.07 



Systematic uncertainties in the double ratios R± and in the CP asymmetries Acp± arise pri- 
marily from uncertainties in signal yields due to the estimate of the peaking backgrounds and from 
the imperfect knowledge of the PDF shapes. The systematic uncertainty associated to peaking 
backgrounds is evaluated by taking into account the uncertainties on their branching fractions and 
by allowing for Poisson fluctuations of their yields in the selected data sample. The estimated yields 
are 29 ± 7 (B-^K~K + K~), 4 ± 4 (B~ ^R-tt+tt-) and O.Olj^ (B~ -» K^tt°K~). Possible CP 
asymmetries up to 30% in their yields are also taken into account. The parameters of the PDFs 
that are fixed in the nominal fit are varied by ±lcr and the difference in the signal yields is taken 
as a systematic uncertainty. 
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The uncertainties in the branching fractions of the channels contributing to the BB back- 
ground have been taken into account. The correlations between the different sources of systematic 
errors, when non-negligible, are considered. An upper limit on intrinsic detector charge bias due 
to acceptance, tracking, and particle identification efficiency has been obtained from the measured 
asymmetries in the processes B~^D°[^K~ir + ]h~~ and B~ — >D(jp±ir~, where CP violation is ex- 
pected to be negligible. This limit (±0.04) has been added in quadrature to the total systematic 
uncertainty on the CP asymmetry. 

5 SUMMARY 

In conclusion, we have reconstructed B~—>D°K~ decays with D mesons decaying to non-CP 
(K~7r + ), CP-even (K~ K + ,ir~ir + ) and CP-odd (K^ir ) eigenstates. We have measured the CP 
asymmetries A C p+ = 0.40 ± 0.15(stat) ± 0.08(syst), A C p- = 0.21 ± 0.17(stat) ± 0.07(syst), and the 
double ratio of branching fractions R + = 0.87 ± 0.14(stat) ± 0.06(syst), P_ = 0.80 ± 0.14(stat) ± 
0.08(syst). These results improve the previous existing measurements from BABAR [Jj. All results 
presented in this document are preliminary. 
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Figure 1: AE distributions of B~^D°h~ candidates, where a charged kaon mass hypothesis is 
assumed for h. Events are enhanced in B — > D°K~ purity by requiring the Cherenkov angle of 
the track h to be within 2a of the kaon hypothesis. Top: B~ — > D°[K~ir + ]K~ ; middle: -B~ — > 
D^p + [i^ _ A' + , 7r~7r + ]i'r~; bottom: B~ — > DQ P _[KgTT°]K~ . Solid curves represent projections of 
the maximum likelihood fit; dashed-dotted, dotted and dashed curves represent the B — > D°K, 
B —* Z?°7r and background contributions. 
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